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In the present letter, with the help of a semi-analytical code we show how taking into account baryonic 
physics in the ACDM model can solve the discrepancy between the numerical predictions of dark matter haloes 
in the ACDM framework and observations from dwarf galaxy scale to clusters of galaxies without the need of 
different forms of dark matter as recently advocated by [1]. Combining well established results, the paper 
shows, for the first time, how the flat profiles of galaxy clusters and correlations between several of their 
main properties, are naturally obtained when baryon physics is correctly taken into account; how the so 
called ’diversity problem” can be solved and how the challenging, extremely low rising rotation curve of 
1C2574 can very well be reproduced when baryon physics is taken into account. We therefore suggest that 
before introducing new exotic features in the standard cosmological model, albeit legitimate, baryonic 
physics should be treated properly to reach agreement with observations. 
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Introduction The ACDM model, while very successful 
[2, 3], present some issues [e.g. 4]. Particularly trouble- 
some is the discrepancy between the flat density profiles of 
dark matter (DM) dominated dwarf galaxies, Irregulars 
and Low Surface Brightness galaxies (hereafter LSBs), 
high surface brightness spiral galaxies in some cluster of 
galaxies, and the cuspy profile predicted by dissipation- 
less N-body simulations [e.g. 5], dubbed cusp/core prob- 
lem [6, 7]. Better understood in terms of the excess of DM 
in the inner parts of the galaxies rather than of the in- 
ner slope, it connects to the Too-Big-To-Fail (TBTF) prob- 
lem [8, 9]. We also mention the large diversity (hereafter 
dubbed “diversity problem”) in the dwarf galaxies rota- 
tion curves (RCs) at odds with hydrodynamic simulations. 


A possible solution to these problems is to assume that the 
dark matter (DM) component is not cold and this leads to a 
wealth of different models (e.g. self-interacting DM model 
[10], SIDM) recently used by [1] to propose a unified solution 
(at all scales) to the small scale problems of the ACDM model 
(the “deficit problem in halos” in their words). 


The study of [1] claims the difficulty or impossibility for 
the ACDM model to explain the rotation curve of 1C2574 
and several other issues, and proposes the SIDM model as 
a possible alternative to the ACDM model. 


Although conceptually interesting, this approach 
presents serious limitations: a) it discards a-priori the 
possibility of solving those problems within the standard 
ACDM model; b) the SIDM’s solutions presented have the 
additional complication of a cross section depending on 


the object mass [11, 12], avoidable by assuming a velocity 
dependent cross section; c) as shown by [13], the diversity 
problem requires additional baryon physics to the SIDM. 

It is thus of fundamental importance to verify whether the 
ACDM model can solve, at all scales, the problems discussed 
in [1], as the SIDM appears unsatisfactory. 

In the present letter we want to address the question: does 
a unified solution exists to the “deficit problem in halos” in 
the ACDM model without invoking a different physics? 


We will closely follow [1] in showing how the “deficit prob- 
lem in halos” is solved when baryonic physics is properly 
taken into account in the ACDM model. 

The plan of this work is as follows: we first present 
our semi-analytical model. Differently from N-body and 
hydrodynamical simulations, the model makes the dif- 
ferent physical contributions clearer and more easily 
disentangled. We then show how it reproduces the clusters 
(e.g, A2537) presented in [14, 15], explains the RCs of some 
peculiar galaxies (e.g., IC 2574, this last fitted by [13] using 
SIDM and baryon physics) and solves the diversity problem 
discussed by [16]. Finally we conclude discussing our results 
comparing them to those of [1]. 


Model: Here we summarize the model described in 
[17-19]. An improvement to the spherical infall models 
(SIM) [20-22], it includes ordered [23] and random angular 
momentum [21, 22], adiabatic contraction [24, 25], the effect 
of dynamical friction of gas and stellar clumps on the DM halo 
[26-30], gas cooling, star formation, photoionization, super- 


nova, and AGN (Active Galactic Nucleus) feedback [3 1, 32]. 

The model considers an initial proto-structure in the lin- 
ear phase containing DM and baryons (diffuse gas) whose 
abundance is given by the “universal baryon fraction” 
fə = 0.167 + 0.01 [3]. Initial conditions are set as in 
[17]. This perturbation expands with the Hubble flow to 
a maximum radius (turn-around) and then re-collapses, 
first in the dark matter component which creates the po- 
tential well baryons will fall in. The equations of motion of 
DM and baryons depend on the angular momentum (AM) 
acquired through tidal torques (ordered AM, L) [23, 33- 
35], on the random AM, j, produced by random velocities 
generated in the perturbation collapse [21], on the dynam- 
ical friction (DF) and on the cosmological constant. The 
calculation of the ordered AM follows the standard Tidal- 
Torque-Theory (TTT), as shown in [17] while the random 
AM, as described in [17], is assigned at turnaround. The 
force between DM and baryonic clumps due to dynami- 
cal friction is evaluated as in [36] and represents an ex- 
tension of [37]. Clumps dimensions, lifetime and charac- 
teristics are described in [19]. For spiral galaxies, they 
are obtained using Toomre’s criterion [see 19]. Their mass 
typically reaches 1% of the disk mass. Clumps density and 
rotation velocity are similar to those found by [38, see their 
Figs. 15 and 16]. 

DM is adiabatically contracted during the collapse 
phase and a more cuspy profile is formed [24, 25]. Assum- 
ing proportionality between DM and baryons initial den- 
sity profiles [39-41, i.e., choosing an NFW profile], Mam 
is obtained by solving iteratively the angular momentum 
conservation equation [23, 24]. 

Baryons cool down due to radiative processes and form 
clumps which collapse to the centre of the halo because of 
DF between baryons and DM, while forming stars. These 
clumps then transfer energy and AM to DM [26, 27], 
increase their random motion and produce a predomi- 
nant outward motion for DM particles, reducing the cen- 
tral density; the density cusp is heated and a core forms. 
This provides the main mechanism of core formation for 
dwarf spheroidals and spirals while for spiral galaxies, this 
mechanism is amplified by the acquired ordered and ran- 
dom AM. 

A similar mechanism contributing to the core forma- 
tion is given by SN feedback [42], in which repeated SN 
explosions flatten the profile. They are compared in de- 
tail in [18]. Gas cooling is treated as a classical cooling 
flow [see 32]; star formation, reionization and SN feed- 
back are included as in [31, 32]. The star formation rate 
W% = asrMst/tayn depends on the star formation effi- 
ciency agp = 0.03 [32], on the dynamical time tayn and 
on the gas mass above a given density threshold Mr, fixed 
to n > 9.3 cm~® as in [43]. The initial mass function 
(IMF) is chosen as Chabrier [44]. At each time step At, 
AM, = WAt stars are generated [see 32, for more details]. 

The specific implementation for SN feedback is given in 
[45] and develop as follows: the energy injected in the in- 


terstellar medium (ISM) by a SN explosion, A Esn is pro- 
portional [see 32] to the typical energy released in the ex- 
plosion, Esn = 10°! erg, to the the number of SN per solar 
mass, nsn = 8 x 1073 M5" (for a Chabrier IMF), to the 
efficiency energy is able to reheat the disk gas chalo = 0.35 
[32] and to AM,, the mass converted in stars. The gas 
reheating produced by energy injection is proportional to 
the stars formed A M,eheat = €aisk AM, where egisk = 3.5 
[32]. The hot gas ejected is proportional to A Esn — AF ot, 
where AF},¢ = 0.5AMyreheattsn Esn is the thermal en- 
ergy change produced by the reheated gas [see 32]. Reion- 
ization reduces the baryon content and takes place in the 
redshift range 11.5-15. 

In the following stage, SN explosions repeatedly eject gas, 
and thus lower stellar density. Feedback destroys the small- 
est clumps soon after a small part of their mass is transformed 
into stars [30]. 

For masses M ~ 6 x 10!4Mo, AGN quenching must be 
taken into account [46]. AGN feedback follows the prescrip- 
tion of [47, 48]. A Super-Massive-Black-Hole (SMBH) is 
created when the star density exceeds 2.4 10° Mo / kpc”, the 
gas density reaches 10 times the stellar density, and the 3D 
velocity dispersion exceeds 100 km/s. Each initial (seed) 
black hole mass starts at 10° Mo. Mass accretion into the 
SMBH and AGN feedback were implemented modifying 
the model by [49] as in [47]. 

Our model does not contain data fitted parameters. Its 
parameters are the same found in simulations (e.g., star 
formation rate, density treshold of star formation, and 
other parameters previously discussed and also found in 
the references given). 

The model demonstrated its robustness predicting be- 
fore simulations the density flattening and shape pro- 
duced by heating of DM, for galaxies [17, 50, 51] in agree- 
ment with subsequent SPH simulations [e.g. 52, 53] [see 
also Fig. 4 of 54, for a direct comparison], and simi- 
larly for clusters of galaxies [55] in agreement with hydro- 
simulations of [56], the inner density slope dependence on 
the halo mass and on the total baryonic content to the to- 
tal mass ratio [18, 50], in agreement with [43]. The slope 
was shown by [18, 50] to depend also on the angular mo- 
mentum. In [19], the stellar and baryonic Tully-Fisher, 
Faber-Jackson and the stellar mass vs. halo mass (SMH) 
relations were shown in agreement with simulations. 

Finally, the correct DM profile inner slope dependence 
on halo mass was explained over 6 order of magnitudes in 
halo mass, from dwarves to clusters[17, 50, 51, 55], a range 
no other model achieved. 

Clusters: The cusp-core problem extends to the cluster of 
galaxies scales: combining weak and strong lensing and stel- 
lar kinematics, the total inner density profile was shown as 
well described by dissipationless N-body simulations at radii 
> 5 — 10 kpc, while DM profiles are flatter than those ob- 
tained in the simulations [14, 15], within a radius of ~ 30 
kpc, typical of the Brightest Cluster Galaxy (BCG) radius. 
The DM profile is characterized by a variation of the slope, 
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FIG. 1. Left panel: Density profile of the total and DM mass for the cluster A2537. The bottom blue (upper green) band represents the DM 
(total mass) density profile determined by [14, 15]. The band in black dotted lines is the DM density profile obtained in this paper. The band 
widths represent the 1 — o uncertainty. The bottom arrow, is the three-dimensional half-light radius of the BCG. The segment with slope 
r—'-13 spans the radial range r = [0.003 — 0.03]r200. Middle panel: inner slope of DM haloes, a ys. the effective radius R. of the BCG 
of the indicated clusters obtained fitting the clusters density profile with the gNFW model. The larger (smaller) error-bars correspond 
to the results by [15] (our model). The dotted lines are the least-square fits. Right panel: core radii, rcore, obtained fitting the clusters 


density profile with a cNFW model vs. Re. 


a = —dlog ppm /dlogr, from cluster to cluster which corre- 
lates with the BCG properties. 

The total and the DM density profiles of MS2137, A963, 
A383, A611, A2537, A1667 and A2390, were determined in 
the aforementioned works. In [14, 15] improved data allowed 
the determination of the stellar mass scale, allowing to pro- 
duce a more physically consistent analysis, reducing the de- 
generacies among stellar and dark mass, and taking into ac- 
count the BCGs homogeneity. 

Two different approaches are available to compare the 
result of our model to the density profiles in [14]. While the 
density profile depends on the virial halo mass, Myir, on the 
baryon fraction, fe = Mp/Mvir, and on the random AM, as 
shown in [55], one can adjust the value of j, so that ppm 
reproduces the observed clusters profiles. Here we rather pre- 
fer the possibility to “simulate” the formation and evolution 
of clusters with similar characteristics to those of [14]: final 
halo mass', Mgco, core radius, etc. 2 

To obtain reasonable agreement between ours’ and ob- 
served clusters, we run simulations till the final halo and 
baryonic masses differ by at most 10% from the observed 
clusters. To get this results the clusters were re-simulated 
50 times. 

Of all the seven cited clusters, in Fig. 1 (left panel) we 
plot the spherically-averaged density profile of A2537 (stud- 
ied by [1]) for the DM halo, BCG stars and their sum (total 
mass) and we refer to the caption for details. The plot shows 
a good agreement between the observations and the model: an 
inner DM profile with almost flat slopes in the case of A2537. 
However, in all the sample the inner slope has an average of 


1 We use Magoo as in [14] 

2 By “simulate” we mean that, as in hydrodynamic simulations, we fix 
the initial conditions and follow the evolution of galaxies, or clusters, 
from the linear to the non-linear phase, until the object formation, 
and its subsequent evolution due to the physical effects previously de- 
scribed. 


~ 0.54, flatter than the NFW profile; a total mass profile close 
to a NFW profile and baryons dominate the profile in the inner 
~ 10 kpe. 

In the middle and right panels we show some correla- 
tions found by [15]. The points with error-bars in the mid- 
dle panel show the value of the inner slope a vs. the effec- 
tive radius, Re, of the BCG for the clusters indicated, ob- 
tained fitting their profiles with the generalized Navarro- 
Frenk-White profile (gNFW). The larger error bars repre- 
sents the result by [15], the smaller ones, our results. The 
right panel represents the core radii, core, vs. Re for the 
same clusters, obtained fitting the density profiles with a 
cored NFW model (CNFW) [15]. Dotted lines are the least- 
square fits. We point out that while the correlations ob- 
served by [15] are re-obtained in our model, simulations, 
usually reaching mass scales of 10'' — 10!? Mo, do not. 
Even the simulations by [47, 56] observe a flattening of the 
inner profile but do not study the correlations. 

Dwarf and LSBs galaxies: We use our model to simulate 
100 galaxies in a ACDM cosmology with similar character- 
istics to the SPARC sample [57], a collection of nearby galax- 
ies high-quality RCs. The stellar mass of the simulated sample 
is in the range M, = 6 x 10°—10!! Mo. Of the galaxies used 
by [1], IC 2574, NGC 2366, DDO 154, UGC 4325, F563-V2, 
F563-1, F568-3, UGC 5750, F583-4, F583-1 are included in 
our larger sample. The observed RCs were compared to the 
most similar simulated galaxies (e.g., with same halo and 
baryonic mass). As an example, in Fig. 2 we show the RC of 
1C2574, the same reproduced by SIDM in [1]. 


Similarly to the case of clusters, we run several simu- 
lations for this galaxy until the final halo and baryonic 
masses differ by at most 10% from the observational data. 
The plot shows the observational RC (dots with error-bars), 
the contribution to the RC given by gas (dashed red line), stars 
(dotted green line), and the total baryonic mass (continuous 
blue line). 


As can be seen here, the RC of [C2574 (advocated to be 
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FIG. 2. Left panel: Rotation curve of IC2574, dot with error bars, from the SPARC catalogue. The solid grey line is the RC of our model. 
The stars contribution to the RC is represented by the dotted green line, the gas disk by the dashed red line, and the total baryonic mass by 
the continuous blue line. This RC is presented for comparison with that shown by [1]. Middle panel: Effects of baryonic physics on the 
relation V2kpe-Vmax. The dashed line represents the expectation if the haloes were all described by a NFW profile. The dots come from the 
galaxies of [16], while the open squares represent the prediction of our model. Right panel: Inner slope of the DM halo vs. /,./Mypaio- The 
dashed line is the result of [43], the dot-dashed and the solid lines are the present paper and the results by [58]. The shaded blue and 
gray regions represent the 1-c scatter in our result and [58], respectively. 


problematic for the scenario of core formation [16]) is very 
well described by the simulated galaxy, as well as the bary- 
onic mass, as seen comparing with [59] or with the SPARC 
mass models [60]. We want to stress that the differences 
between the [C2574 RC’s baryon contribution (gas, stars) 
in [1] from that of SPARC or [59] arise from the differ- 
ent DM properties. The halo mass of the simulated galaxy 
host is 1.8 x 101? Mo, as obtained in the SPARC sample 
or from [60]. 

Recently, [61], showed how the RC of IC2574 can be 
naturally obtained taking into account SN feedback, with 
a similar approach to ours. 

Diversity: Despite the fact that the RCs of dwarf galaxies 
are on average cored, individual fits to galaxy RCs show inner 
slopes ranging from a ~ 0 to cusps, while for cored profiles, 
the central densities can differ by a factor of 10 for galax- 
ies inhabiting similar halos [62] and the situation becomes 
more complicated at higher masses. For many objects [63] 
found cored and cuspy profiles in dwarves which are sim- 
ilar while [64] a tendency to flatter profiles in less massive 
galaxies. 

Such diversity was quantified by [16] comparing the circu- 
lar velocity at 2 kpc V2kpc, with a fixed value of the maximum 
of circular velocity (Vmax). For 50 < Vmax < 250 km/s there 
is a scatter of a few in V2kpc- [13] studied the problem in the 
SIDM scenario, and found that SIDM alone cannot explain 
the scatter, since the resolution requires baryonic physics 
must also be taken into account. 

In Fig. 2, following [13] and [16], we plot (left panel) 
Vakpe versus the maximum of the circular velocity Vmax- 
The dashed line is the expectation for a NFW density pro- 
file. The dots are the observed values as in [16] and the open 
squares are the prediction of our model. The sample used 
by [16] contains galaxies with stellar mass in the range 
Mstar ~ 107 — 10!! Mo. In order to compare their re- 


sults with our model, we chose simulated galaxies with 
stellar mass in the range of [16] and we calculated Vzķpc 
and Vmax. 

As the dots show, at fixed Vmax, the scatter in V2kpc can be 
as large as a factor of four. Such scatter cannot be explained 
by the ACDM model, as it produces cuspy and self-similar 
halo density profiles, with a single parameter (concentration 
parameter or halo mass), in contrast to the cores displayed 
by many dwarves. In the ACDM model, the much larger 
amount of DM in the halo cusp than the baryons “freezes” 
the scatter in V2kpc, produced, conversely, by the spread in 
the baryon distribution. Baryon physics heats DM and enlarge 
the galaxies, reducing the inner DM content. As shown in 
[19, 50], the inner slope of the halo density profile is mass de- 
pendent. Milky Way sized galaxies tend to have cuspy profiles 
while dwarf sized galaxies cored profiles. Ultra faint dwarf 
galaxies tend to be more cuspy than dwarves. Therefore, the 
scatter seen in Fig. 2 originates from the mass dependence of 
the core formation process and the effects of environment, as 
described in [51]. 

Our model successfully recovers the scatter and distri- 
bution of the RCs shapes because baryon physics gives rise 
to different responses in the halo of simulated galaxies. 
The right panel represents the inner slope of the DM halo 
obtained in [43], dashed line, in the present paper, dot- 
dashed line, and in [58], solid line. All the curves were ob- 
tained as in [43] by fitting the DM profile with a power law 
in the radial range 0.01 < r/Ryir < 0.02, being R,;i, the 
virial radius. The shaded blue and gray region represents 
the 1-o scatter in our result and [58]. 

The plot shows that the core formation mechanism and 
a are strongly dependent on M, /Mhalo ° with a minimum 


3 The correlation between a and Mx /Mhalo can be expressed in terms of 


value of a at masses M,./Mpaio ~ 107? corresponding to 
M, ~ 108 Mo [18, 43, 58] due to the maximum effects of 
baryon physics. For smaller masses the profile steepens 
because of the relative decrease of stars (ratio M, /Myaio)- 
Since the profile tends to steepen for M, < 10°M. and 
Vmax is proportional to V., we should expect a self-similar 
behaviour, similar to the NFW RCs, as observed in Fig. ??. 
For M, > 10°Mo, the increase in stellar masses gives rise 
to a deepening of the potential well and a reduction of the 
effects of baryon physics, with a consequent steepening of 
the profile. 

Conclusions: The ACDM model exhibits some problems 
at small scales, and in particular predicts an excess of DM 
in the central parts of galaxies and clusters. In this letter, 
we showed that a unified solution to the problem can be ob- 
tained within the ACDM framework without introducing dif- 
ferent forms of DM [as done in |, instead]. With a semi- 
analytic model, we simulated the clusters studied by [15] and 
compared the density profiles with those they obtained :those 
profiles were re-obtained correctly by our model. We dis- 
played one of those density profiles (Fig. 1). We then sim- 
ulated a sample of galaxies similar to the SPARC compila- 
tion, containing the galaxies studied by [1], finding again a 
good agreement with data, as shown in the case of one of 
the most complicated galaxy RCs to reproduce, namely that 
of 1C2574 (Fig. 2). Finally, we studied the “diversity” prob- 
lem using the simulated galaxies, and comparing their Vaiyc, 
for given values of Vmax With the compilation in [16]. We 
show that baryon physics gives rise to RCs very different from 
each other, due to the dependence of the RC from their total 
and stellar mass, together with environment. This explains the 
scatter in the Voxpe-Vmax plane. 
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